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In the search for the mechanism of high-temperature superconductivity, intense research has 
been focused on the evolution of the spin excitation spectrum upon doping from the antifer- 
romagnetic insulating to the superconducting states of the cuprates.^"^ Because of technical 
limitations, the experimental investigation of doped cuprates has been largely focused on low- 
energy excitations in a small range of momentum space.^"^^ Here we use resonant inelastic 
x-ray scattering^^"^^ to show that a large family of superconductors, encompassing underdoped 
YBa2Cu408 and overdoped YBa2Cu307, exhibits damped spin excitations (paramagnons) with 
dispersions and spectral weights closely similar to those of magnons in undoped cuprates. The 
comprehensive experimental description of this surprisingly simple spectrum permits quantita- 
tive tests of magnetic Cooper pairing models. A numerical solution of the Eliashberg equations 
for the magnetic spectrum of YBa2Cu307 reproduces its superconducting transition temperature 
within a factor of two, a level of agreement comparable to Eliashberg theories of conventional 
superconductors.^^ 



The twenty-fifth anniversary of the discovery of high-temperature superconductivity is 
approaching without a clear and compelling theory of the mechanism underl5dng this phe- 
nomenon. After the discovery of an unconventional (d-wave) sjonmetry of the Cooper pair 
wave function in the copper oxides, the thrust of research has been focused on the role of 
repulsive Coulomb interactions between conduction electrons, which naturally explain this 
pairing S5niimetry. However, since even simple models based on repulsive interactions have 
thus far defied a full solution, the question of whether such interactions alone can generate 
high-temperature superconductivity is still open. A complementary, more empirical approach 
has asked whether antiferromagnetic spin fluctuations, which are a generic consequence of 
Coulomb interactions, can mediate Cooper pairing in analogy to the phonon-mediated pairing 
mechanism in conventional superconductors.^ This scenario requires the existence of well- 
defined antiferromagnetic spin fluctuations in the superconducting range of the cuprate phase 
diagram (for mobile hole concentrations 5% < p < 25% per copper atom), well outside the 
narrow stability range of antiferromagnetic long-range order (0 < p < 2%). 

An extensive series of experiments using inelastic spin-flip scattering of neutrons has indeed 
revealed low-energy spin fluctuations in doped cuprates.^ Signatures of coupling between spin 
and charge excitations have also been identified,^ and evidence has been reported that this 
coupling is strong enough to mediate superconductivity in underdoped cuprates."^ However, 
for cuprate compounds hosting the most robust superconducting states, namely those that 
are optimally doped to exhibit > 90 K, inelastic neutron scattering (INS) experiments 
have thus far mostly revealed spin excitations over a narrow range of excitation energies 
£ ~ 30 — 70 mey wave vectors Q covering only ~ 10 % of the Brillouin zone area around 
the antiferromagnetic ordering wave vector Q^, and temperatures T < T^.^~^^ The energy- 
and momentum-integrated intensity of these excitations constitutes only a few percent of 
the spectral weight of spin waves in antiferromagnetically ordered cuprates,^'^® and is thus 
clearly insufficient to support high-T^. superconductivity. Although the recent discovery of a 
weakly dispersive magnetic excitation in the model system HgBa2Cu04+5 with £ ~ 50 meV may 
account for some of the missing spectral weight, the apparent weakness of antiferromagnetic 
fluctuations in optimally doped compounds has been used as a central argument against 
magnetically mediated pairing scenarios for the cuprates.^° 

This picture is, however, strongly influenced by technical limitations of the INS method that 
arise from the small cross section of magnetic neutron scattering in combination with the weak 
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primary flux of currently available high-energy neutron beams. Because of intensity constraints, 
even the detection of undamped spin waves in antiferromagnetically ordered cuprates over 
their full band width of ~ 300 meV has required single-crystal samples with volumes of order 
10 cm"^, which are very difficult to obtain. ^^"^'^ Doping further reduces the intensity of the INS 
profiles and exacerbates these difficulties. 

Here we take advantage of recent progress in the development of high-resolution res- 
onant inelastic x-ray scattering (RIXS) in order to explore magnetic excitations in a wide 
energy-momentum window that has been largely hidden from view by INS. Experiments on 
undoped cuprates have shown that RIXS with photon energies at the Cu L3 absorption edge is 
sensitive to single-magnon excitations ^"^"^^ by virtue of the strong spin-orbit coupling of the 
2p3/2 core-hole intermediate state,^^'^^ in excellent agreement with INS measurements.^"^'^^ 
Initial RIXS experiments on doped La2_;cSrj^Cu04 have revealed a complex two-component 
lineshape with excitations extending up to about 350 meV; the presence of two branches was 
attributed to phase separation between the superconducting state and a competing state with 
incommensurate spin and charge orde r.^^ We have applied the same method to YBa2Cu306+^ 
(YBCOe+x), Ndi 2Bai 8CU3O6+X (NdBCOe+x), and YBa2Cu408, which are much less affected by 
doping-induced disorder and do not show phase separation.^^ 

Although one cannot reach with RIXS, one of its great advantages over INS is that it 
allows measurements of magnetic excitations with sizable intensity over much of the accessible 
reciprocal space, even on very small sample volumes. The results presented here have been 
obtained on thin films and on millimeter-sized single crystals far below the volume requirements 
of INS (see Supplementary Information). Figure 1 shows a sketch of the scattering geometry of 
our experiment, as well as RIXS spectra obtained on undoped antiferromagnetic NdBCOe (Fig. 
1 c,d) and underdoped superconducting NdBCOy (Fig. 1 f,g), for incident photon polarizations 
in and out of the scattering plane (tt and a geometries) and a momentum transfer q// along 
the reciprocal-space a* direction corresponding to 0.37 reciprocal lattice units (r.l.u.). For both 
scattering geometries, the spectra exhibit an intense peak located around 1.7 eV energy loss 
that arises from optically forbidden dd excitations (that is, transitions of the unpaired hole of 
Cu^"*" from the dx2_^2 to other d orbitals)^^ on top of a continuum of charge-transfer excitations. 
At lower energies, we can see in both compounds an inelastic feature centered around 250 
meV and a resolution-limited elastic peak. 

For undoped antiferromagnetic NdBCOg, we can decompose the response in the mid- 
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infrared (MIR) region of the spectra (Fig. le) for energy losses below ~ 500 meV following 
the method employed in Ref. 27. In the n scattering geometry, this leads to (i) an intense 
resolution-limited peak around 250 mey (ii) a high-energy tail of this peak centered around 
400 mey and (iii) a weak low-energy contribution around 100 meV In the cr scattering 
geometry, feature (i) is strongly suppressed. This polarization dependence allows us to assign 
this feature to a single-magnon excitation, in agreement with theoretical considerations^^'^^ 
and previous investigations on other cuprates.^'*"^^'^'^ This assignment is further confirmed by 
the disappearance of the MIR inelastic response as the incident photon energy is moved away 
from the Cu-Lg ab sorption edge (not shown). The weak features (ii) and (iii) are associated 
with higher-order spin excitations and lattice vibrations (single and multiple phonon excitations 
that are not individually resolved), respectively. 

The energy of the single-magnon feature in NdBCOg depends strongly on q// (left panel 
of Fig. 2). Since the lineshape of the MIR spectrum is q// -independent, the fitting procedure 
described above can be used to obtain the magnon dispersion relation. The result displayed in 
Fig. 3a is clearly different from the sinusoidal magnon dispersion observed in La2Cu04.^^'^^ 
Indeed, due to the presence of two antiferromagnetically coupled CUO2 layers per unit cell, 
two magnon branches with different spin precession patterns are expected^^"'^'' for NdBCOg: an 
acoustic branch with intensity oc sin^(^), where is the momentum transfer perpendicular 
to the CUO2 layers and d is the spacing between the two CuOa layers, and an optical branch 
with intensity oc cos^(^). The result of the calculation of the relative intensities of these two 
branches for our scattering geometry is displayed in the inset of Fig. 3a. It shows that close to 
the Brillouin zone center the gapped optical branch dominates. Fitting our data using the spin- 
wave dispersion calculated for a bilayer in the framework of a simple Heisenberg model,^^"^'^ 
we obtain J//=133 i 2 meV and J^=12 ± 3 meV for the intra- and inter-layer exchange 
constants, respectively, in excellent agreement with the values obtained in antiferromagnetic 

We now turn to the doped systems. Figure 2 provides a S5niopsis of the experimental 
spectra for all systems investigated here, which span a wide range of doping levels: undoped 
NdBCOg, strongly underdoped NdBCOy and YBCO5 g (T, = 65 K and 61 K, respectively), weakly 
underdoped YBa2Cu408 (T, = 80 K), and weakly overdoped YBCO7 (T, = 90 K). Because of 
their stoichiometric composition and electronic homogeneity, the latter two compounds have 
served as model compounds in the experimental literature on high-T^ superconductivity, but 
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apart from the "resonant mode" that appears in YBCO7 below T^^ no information has been 
available on their magnetic excitation spectra. Fig. 2 shows that broad MIR features are present 
at all doping levels in the same energy range as the single-magnon peak in undoped NdBCOg. 
Since these features obey the same polarization dependence as the magnon mode (Figs. 1 f,h), 
they can be assigned to magn etic excitations. 

Based on the metallic nature of the doped cuprates, one expects strong damping of magnetic 
excitations in the Stoner continuum of incoherent electron-hole excitations. We have therefore 
fitted these spectra to Voigt profiles that are the result of the convolution of the Lorentzian 
lineshape of excitations with finite lifetime with the Gaussian resolution function. Since the 
fits yield excellent agreement with this simple profile (solid lines in Figs. 1 and 2), we can 
accurately extract the energies and half-widths-at-half-maximum (HWHM) of the magnetic 
excitation as a function of the transferred momentum (Fig. 3b, c). Before discussing these results, 
it is interesting to compare them to those previously reported on underdoped La2_;<:SrxCu04.^'^ 
As no double-peak structure is observed here, our data confirm the absence of phase separation 
in the YBCO and NdBCO families of compounds.^^ 

The magnetic excitation energies of NdBCOy, YBCOe.e, and YBCO7 at the Brillouin zone 
boundary are nearly identical to the ones found in NdBCO^. This implies that the in-plane 
exchange constant J// is not as strongly renormalized with doping as previously suggested 
based on extrapolation of lower-energy INS data,^°"^^ in agreement with the most recent INS 
findings on YBCO5 5 using high energy neutrons.^^ Upon approaching the T point (q// = 0), we 
observe that the dispersion in the doped compounds is steeper than in NdBCOg. Despite the 
fact that r and are no longer equivalent in the absence of magnetic long-range order, the 
RIXS data obtained on YBCO5 5 nicely extrapolate to the low-energy "hour glass" dispersion 
around previously extracted from INS data on samples prepared in an identical manner.^"* 
Since the slope of the upward-dispersing branch close to the magnetic zone center is flatter 
than the one of the antiferromagnetic spin-wave dispersion,^^'^'^ this indicates the presence of 
an inflexion point in the dispersion of magnetic excitations in the doped compounds. We note 
that the situation seems somewhat different for YBa2Cu408, where the energy of the magnetic 
excitations close to the zone boundary is significantly reduced compared to the other systems 
(~ 210 meV instead of ~ 300 meV), and the dispersion is flatter than in NdBCOs- This may 
reflect different values of J// and in this system, and deserves further investigation. The 
intrinsic HWHM of ~ 200 meV of the inelastic signal extracted from our data (Fig. 3c) is much 
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larger than the instrumental resolution and comparable to the magnon energies, indicating 
strong damping by Stoner excitations. The damping rate does not change substantially with 
qil and with doping. Finally, while it is not possible yet to obtain absolute magnetic intensities 
from RIXS data, we can extract the relative intensity of the RIXS profiles by integrating the 
inelastic signal in the MIR region (see Supplementary Information). Remarkably, the integrated 
intensity obtained in this way is conserved upon doping from the antiferromagnetic insulator 
to the slightly overdoped superconductor. 

We have thus demonstrated the existence of paramagnons, Le^ damped but well-defined, 
dispersive magnetic excitations, deep in the Stoner continuum of cuprates with doping levels 
beyond optimal doping. Their spectral weights are similar to those of spin waves in the 
undoped, antiferromagnetically ordered parent material. 

In order to obtain insight into the origin of this surprising observation, we have performed 
exact-diagonalization calculations of the t — J Hamiltonian with 3 jt = 0.3 on finite-sized 
clusters, following the method proposed in Refs. 35-37. We used clusters with 18 and 20 spins 
arranged in various shapes to map the (H,K,0) reciprocal plane with sufficient momentum 
resolution (see Supplementary Information). The resulting spectra were convoluted by a 
Gaussian function with HWHM 0.1 1 in order to simulate the experimental resolution function. 
The calculated imaginary part of the spin susceptibility is shown in Fig. 4a for different hole 
concentrations. In the undoped case, one can observe two peaks in the imaginary part of 
the spin susceptibility: an intense peak corresponding to the single-magnon excitation, and a 
weaker feature at higher energy corresponding to higher-order processes. The single-magnon 
peak clearly disperses, although its energy is slightly above the value expected from linear 
spin-wave theory (dotted lines), due to the finite size of the cluster. In the inset of Fig. 4a, we 
present the energy-integrated magnetic intensity in the (H,K,0) reciprocal plane. As expected, 
the intensity for in the u ndoped situation is strongly peaked at Q^, and rather uniformly 
distributed in the rest of the plane. As holes are added to the clusters, the magnetic spectral 
weight strongly decreases only around Q^, but remains essentially constant everywhere else. 
This can also be seen in Fig. 4b, where the imaginary part of the energy-integrated spin 
susceptibility is plotted as a function of doping. The experimental and numerical data are thus 
in excellent agreement. 

Armed with essentially complete knowledge of the spin fluctuation spectrum, and motivated 
by the agreement with the numerical data for the t — J model, we now estimate the super- 
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conducting transition temperature generated by a Cooper pairing mechanism in which the 
experimentally detected spin fluctuations play the role of bosonic glue. To this end, we have 
self-consistently solved the Eliashberg equations using the vertex function of the t — J model^® 
and the experimental spin fluctuation spectrum of YBCO7 (inset in Fig. 4c), without adjustable 
parameters. Details are given in the Supplementary Information. Figure 4c shows the outcome 
of the calculation. The resulting of 170 K is similar to the maximum observed in the 
cuprates and to another recent estimate based on the comparison of INS and photoemission 
data on underdoped ¥6005+^."^ The agreement with the experimentally observed of YBCO7 
is satisfa ctory in view of the neglect of vertex corrections beyond Eliashberg theory and 
other simplifying assumptions in the calculation, and in view of a similar level of quantitative 
agreement reported for Eliashberg calculations of conventional low-T^. superconductors.^^ In 
order to resolve the relative contribution of the strongly doping dependent region around 
and the high-energy part of the spin susceptibility revealed by RIXS to the pairing strength, 
we have introduced momentum cutoffs around in the calculation (Fig. 4c). Clearly, both 
low- and high-energy spin fluctuations contribute substantially to the pairing. In underdoped 
cuprates with lower gaps and stronger antiferromagnetic correlations, vertex corrections are 
expected to become essential,^^ reducing the values calculated here. These corrections 
suppress coupling to magnons in the vicinity of Q^, leaving the contribution of the higher-ene 
rgy excitations as a vital source of pairing. 
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Fig. 1: a and b) Schematics of the scattering geometry. Typical RIXS spectra of undoped antiferromag- 
netic Ndi 2Bai gCusOe (c, d, e) and superconducting underdoped Nd^ 2Bai.8Cu307 (f, g, h), obtained at 
T = 15 K for qj 1=0.27 r.l.u., in both n (black squares) and cr (red squares) scattering geometries. 



AF-NdBCO, NdBCO, YBCO,, Y124 YBCO, 

6 7 6.6 7 
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Fig. 2: RIXS response of undoped Ndi sBa^ gCugOe (NdBCOe), underdoped Nd^ sBai gCugOy (NdBCOy), 
YBasCugOg 6 (YBCOe.e), YBa2Cu408, and slightly overdoped YBasCugOy (YBCO7) for various momen- 
tum transfers (see Fig. 1), at T = 15 K. Note that the spectra for NdBCOy were recorded for 5=20°, 
corresponding to q// = 0.18 r.l.u. 
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Fig. 3: a) Experimental magnon dispersion along 100 direction in AF Nd^ 2^31 8CU3O5 at T = 15 K, 
fitted using the spin- wave dispersion of a bilayer from Ref. 21 (thick red line). The dashed lines are the 
acoustical (black) and optical (red) spin wave dispersions calculated using the fitting parameters. The 
grey area represents our energy/momentum resolution. Inset: relative intensity of the acoustical and 
optical magnon for our scattering geometry, b) Experimental magnon dispersion along 100 direction 
in antiferromagnetic Ndi 2Bai gCusOe, underdoped Ndi 2Bai gCusOy, YBa2Cu306 YBa2Cu408 and 
YBa2Cu307 at T = 15 K. Low-frequency INS data recorded along the 100 direction from (tt, ti) for 
YBa2Cu306 5 have been added."'^ Lines are guides to the eye. c) HWHM of magnetic excitations in 
Ndi 2Bai 8CU3O7, YBa2Cu305 5, YBa2Cu408 and YBa2Cu307. d) Integrated inelastic intensities (see also 
Supplementary Information). 
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Fig. 4: a) Imaginary part of the spin susceptibility resulting from exact diagonalization of the t — J 
model with J/ 1 = 0.3 on small clusters. Results from 18-site and two 20-site clusters are combined to 
cover the F —X line (see Supplementary Information) . The spectra are broadened using a Gaussian 
with HWHM= 0.1 1. Dashed lines correspond to the linear spin wave dispersion with the same J as 
used in numerics. Insets: Brillouin-zone-maps of energy-integrated x" with a common color scale, 
b) Energy- integrated x" of the 20-site cluster, normalized to the number of electrons on the cluster, 
as a function of doping. The seven accessible non-equivalent q-vectors for this cluster are shown in 
the inset, c) Superconducting transition temperature resulting from the Eliashberg calculation for the 
experimentally determined spin excitation spectrum of YBCO7 displayed in the inset. The red marks 
indicate the influence of successive momentum-space cutoffs limiting the maximum distance from 
of the q-vectors for which Xqi^) is included in the calculation. 
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Supplementary Information for 
"Intense paramagnon excitations 
in a large family of high-temperature superconductors" 

Sample Preparation and Characterization 

We used untwinned YBa2Cu305 g and YBaaCUgOy single crystals of volume about 3 x 3 x 
0.5 mm'', with superconducting transition temperatures of 61 K and 90 K (AT^ =2 K), 
respectively, determined for each crystal by SQUID magnetometry (see details in Methods 
section of Ref. 1). The YBa2Cu408 single crystals were obtained using KOH flux under ambient 
pressure in a box furnace. The source material was conventionally prepared polycrystalline 
YBaaCUgOy mixed with CuO in a molar ratio of 1:1. The details are described in Refs. 2,3. The 
crystal size was 0.8 x 0.8 x 0.2 mm^, and = 80 K. Nd^ sBa^ 8Cu05+^ films of thickness 100 
nm were deposited on SrTiOs (100) single crystals by diode high-pressure oxygen sputtering. 
The undoped Ndi 2Bai gCuOg film was obtained by annealing as-grown Ndi 2Bai gCuOy samples 
in Argon atmosphere (10 mbar) for 24 hours. Details of the film growth and characterization 
can be found in Ref 4. 

Table I Usts the lattice parameters and transition temperatures of all samples. 

Measurements 

The RIXS measurements were performed at the ADRESS beam line^ of the Swiss Light 
Source (Paul Scherrer Institute, Switzerland) using the high-resolution SAXES spectrometer.^ 
The resonant conditions were achieved by tuning the energy of the incident x-ray to the 
maximum of the Cu L3 absorption peak, around 931 eV The total momentum transfer is 0.855 
A~^, which allows one to cover about 85% of the first Brillouin zone along the [100] direction 
(see Supplementary Figure 2). Momentum transfers are given in units of the reciprocal lattice 
vectors a*, b* and c* where a*=27i/a, b*—27i/b, and c*—2n/c. (See Table I for values of a, 
b, and c for each sample). The total energy resolution was about 130 meX and the exact 
position of the elastic (zero energy loss) line was determined by measuring, for each transferred 
momentum, a non-resonant spectrum of polycrystalline graphite. All data where recorded at 
r = 15 K, each spectrum being the result of 30 to 180 minutes total accumulation (sum of 
individual spectra of 5 min). 
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Data Analysis 

The data on the parent compound were analyzed using resolution-limited Gaussian line- 
shapes, as in Ref. 7. Due to the broadening of the lineshape with doping, this analysis does 
not hold an5miore for the other systems. The RIXS intensity is proportional to the djmamical 
structure factor S(Q,co) = [1 + nico,T)]x'XQ,<o), where [1 + n(co,r)] = (1 - e'^BT)-^ is 
the Bose factor (step function at T = 15 K) and x XQ, ^) is the imaginary part of the spin 
susceptibility. To take into account the finite lifetime of the excitations in the doped systems, 
and causality (that imposes x"iQ.> ^) to be an odd function of energy transfer co), we modeled 
x"(.Q., ^) using an antisymmetrized Lorentzian function, centered on a> = coq, and of HWHM 



X'XQ, CO) = 



(1) 



L (CO - coq)2 + (r^ ) (CO + coq)2 + (r^ ) J ' 

convoluted with the Gaussian resolution function (HWHM = 65 meV). Note that in the case of 
doped systems, we are not able to separate the vibrational contribution to the inelastic response 
from the magnetic one. This may result in a slight, systematic overestimate of the intrinsic 
HWHM of the magnetic excitation plotted in Fig. 3c in the main text. As in the undoped case, 
this contribution does not account for more than 10 % of the total inelastic spectral weight 
(the inelastic intensities plotted in Fig. 3d all include this vibrational contribution). 

In Supplementary Figure 1, the x-ray absorption spectra close to the Cu Lg absorption 
edge (normal incidence, n polarization) are shown for each of the investigated compounds. 
These measurements were taken on beamline ID08 at the European Sjmchrotron Radiation 
Facility (ESRF). Since these absorption profiles are essentially identical between the different 
samples, the self-absorption corrections are not at all critical here. The relative change of the 
intensities due to self-absorption obtained by doing calculations of the type given in Ref. 8 are 
smaller than 3 %. The high-energy dd excitation peak (Fig. 1 b,c in the main text) provide an 
additional calibration standard. 



Model calculations 

Supplementary Figure 3 shows the cluster shapes used in the numerical computation of the 
d5niamical spin susceptibility of the t — J model. The combination of different cluster shapes 
displayed in panel a provides access to the momentum points in panel b. 
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The Eliashberg calculations were performed with the vertex function of the t — J model, 
^Yq + (^fc-q + ^fc)/2 where = — 2t(cosfcjf + cosfc^,) — 4t'cosfc^cosfc^ represents the bare 
electronic dispersion and Yq = (cosq^ + cosqy)/2.^ The values of the model parameters used 
in our calculation are t = 400 meX t' = — 1/3, and J — 130 meV As input we used the 
simple formula Xqi^) = Vq/i^l — co^ — ir^co), which provides an excellent description of 
both the experimental data and the spin susceptibility resulting from the exact-diagonalization 
calculations. The parameters co^ and denote the paramagnon dispersion and damping (full 
width at half maximum), respectively. We take cOq = 2J^(1— Yq)(l + Yq + cOq/SJ^) with a 
spin gap coq at Q^. For underdoped YBCO^+j^, coq can be directly extracted from INS data. For 
YBCO7, where normal-state INS data are not available, the resonant mode below implies 
cOq ~ 60 meV This value is required to reproduce the experimentally observed resonant mode 
energy of 40 meV^'^ following method described in Ref. 11. Further, rjq = A(l — Yq) with A 
chosen such that x" obeys the corresponding sum rule, as required by the RIXS and INS data. 

Supplementary Figure 4 gives an overall picture of as a function of the damping r (chosen 
to be independent of q) and the spin gap at Q^. Apparently, non-damped spin waves with 
a suitably high energy at are optimal as the pairing mediators. decreases when the 
damping increases and/or when the spin gap deviates from its optimum value, coq ~ 100 meV 
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Sample 


a 


b 


c 


Transition Temperature 


NdiaBaigCugOg 


3.905 A 


3.905 A 


11.85 A 


Tjv > 300 K 


Ndi.aBai.gCugOy 


3.885 A 


3.92 A 


11.7 A 


= 65 K 


YBaaCugOg.e 


3.82 A 


3.87 A 


11.7 A 


= 61 K 


YBa2Cu408 


3.84 A 


3.87 A 


27.25 A 


r, = 80 K 


YBasCUgOy 


3.817 A 


3.884A 


II.68IA 


= 90 K 



Table I: Supplementary Material Lattice parameters and transition temperatures of the samples inves- 
tigated here. 
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Fig. 5: Supplementary Material Cu L3 edge X-ray absorption spectra of all the investigated compounds 
measured with total electron yield (normal incidence, n polarization). Spectra have been vertically 
shifted by 0.5 for clarity. 
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Fig. 6: Supplementary Material Sketch of the accessible reciprocal space with RIXS experiments at the 
Cu L3 edge. The black line represents the region explored in this work. In red we have represented the 
typical region of the reciprocal space studied INS in doped YBCO. 
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Fig. 7: Supplementary Material a) Cluster shapes used to compute the spin fluctuation spectra of Fig. 
4 in the main text, b) Non-equivalent q-vectors accessible to calculations on the different clusters. 



23 




10 20 50 100 200 500 1000 

damping r [meV] 



Fig. 8: Supplementary Material Contour plot of the superconducting transition temperature resulting 
from the Eliashberg calculation as a function of the spin gap at wave vector and the damping 
parameter P. The values appropriate for YBCO7 are marked in blue. 
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